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ABSTRACT 

We present 3D hydrodynamical simulations of the superbubble M17, also known as 
the Omega nebula, carried out with the adaptive grid code yguazu-A, which includes 
radiative cooling. The superbubble is modelled considering the winds of 11 individual 
stars from the open cluster inside the nebula (NGC 6618), for which there are estimates 
of the mass loss rates and terminal velocities based on their spectral types. These stars 
are located inside a dense interstellar medium, and they are bounded by two dense 
molecular clouds. We carried out three numerical models of this scenario, considering 
different line of sight positions of the stars (the position in the plane of the sky is 
known, thus fixed). Synthetic thermal X-ray emission maps are calculated from the 
numerical models and compared with ROSAT observations of this astrophysical object. 
Our models reproduce successfully both the observed X-ray morphology and the total 
X-ray luminosity, without taking into account thermal conduction effects. 

Key words: ISM: bubbles — X-rays: ISM — Galaxy: open clusters and associations: 
individual (M17, NGC 6618) — stars: winds, outflows — methods: numerical 



1 INTRODUCTION 

, The mechanical luminosity of a stellar wind {\M m V^; where 
' Mm is the mass loss rate and Voa is the terminal speed of 
the wind) is typically less than 1 per cent of the stellar ra- 
diative luminosity. The fast winds of massive stars have, 
however, a great influence to their surrounding interstellar 
medium (ISM). They sweep up the ISM creating a vari- 
ety of structures, from small bubbles around single stars 
to large superbubbles around OB associations. The interior 
of these bubbles or superbubbles contains shock-heated gas 
(with temperatures in excess of 10 6 K), thus emits strongly 
in X-rays, while the outer shell is cooler and bright in opti- 
cal emission lines (Weaver et al. 1977; McCray & Kafatos, 
1987; Mac Low & McCray, 1988; Chu et al. 1995). 

In the last decades, observations from old X-ray satel- 
lites such as the ROSAT, and newer ones, such as the Chan- 
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dra and XMM-Newton observatories, have found diffuse X- 
ray emission inside OB massive star clusters. The diffuse X- 
ray region is associated to hot gas contained into a superbub- 
ble produced by the interaction of the individual OB stellar 
winds with the surrounding environment. Nevertheless, ex- 
tremely massive stars from the cluster would rapidly evolve 
into supernovae (SNe) explosions, in which case the super- 
bubble will be produced by the combined action of stellar 
winds and supernovae. In order to isolate the X-ray emission 
from stellar winds alone we can study very young superbub- 
bles, where SNe explosions are still not taking place. 

The superbubble M17 contains the massive stellar clus- 
ter NGC 6618, with more than 800 stellar sources (Broos et 
al. 2007). This cluster has 100 stars earlier than B9 (Lada 
et al. 1991) and an average age of ~ 1 Myr (Hanson et al. 
1997). With such a short age it is unlikely to have produced 
any SN explosion; in a cluster with 100 O type stars the 
first SN explosion would occur after ~ 4 Myr (Knodlseder 
et al. 2002). Therefore, M17 provides an ideal laboratory to 
observe and study the diffuse X-ray emission as produced 
by stellar wind collisions alone. This kind of superbubbles 
are called 'quiescent superbubbles'. 
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Figure 1. Image of the omega nebulae from the Digital Sky Sur- 
vey (DSS). It is illuminated by the massive stellar cluster NGC 
6618, the centre of which is located behind substantial obscuring 
material, with A n ^ o mag. 



This superbubble is highly asymmetric due to the inter- 
action with the edge of a massive molecular cloud, M17SW, 
located to the West of M17 (Lada et al. 1974). In this direc- 
tion the bubble encounters resistance, while it can expand 
more freely to the East. This asymmetry is evident in a 
large-scale map of M17 at A 21 cm (Felli et al. 1984), where 
two intersecting clouds or 'bars' are observed in emission. 
These have been called the northern (N) and southern (S) 
bars. The projected size of each bar is of ~ 5.7 pc, forming 
an angle of ~ 45° with each other on the plane of the sky 
(Felli et al. 1984; Brogan & Troland 2001). 

The centre of the stellar cluster NGC 6618 consists of 
a 'ring' of seven O-type stars (Hanson et al. 1997) which 
are located between the two bars. These stars must be the 
dominant sources of ionising photons on the nebula, in fact 
diffuse X-ray emission is observed to be well confined within 
the external shell of the optical nebula (Dunne et al. 2003) . 
Figure[T]shown a DSS image of this object. The ring of seven 
O-type stars is located on the West, where is substantial 
obscuring material, with A v ~ 8 mag (Hanson et al. 1997). 

The stellar winds of these seven stars strongly con- 
tribute to the form ation of a common 'stellar cluster wind' 
|Canto et al.ll200oT ). The collision of these winds with the 
surrounding inhomogeneous ISM produces the diffuse X-ray 
emission observed by ROSAT and Chandra satellites. 

In ROSAT images, M17 is observed as a blister-like 
structure, with an overall diameter of ~ 10 — 12 pc, if a 
distance of 1.6 kpc is considered (Nielbock et al. 2001). 
The emission peak is located inside the stellar cluster. This 
suggested that the cluster is the origin of the diffuse X- 
ray emission, and was later confirmed by Chandra obser- 
vations (Townsley et al. 2003), where the angular resolution 
allowed to see the individual sources. M17 was observed with 
ROSAT Position Sensitive Proportional counter (PSPC) for 
6.7 ks in the 0.1 - 2.4 keV band (Dunne et al. 2003). They 
found a total diffuse X-ray luminosity of ~ 2.5 x 10 33 erg s _1 
with kT = 0.66-0.78 keV (T = 7.7 - 9.1 x 10 6 K) and 
N H = (1-5) X 10 21 cm" 2 . 

Several models have tried to reproduce the total X-ray 
luminosity Lx from superbubbles, often taking into account 



the effects of the electron thermal conduction. However, the 
Lx reproduced by models where thermal conduction is in- 
cluded (Weaver et al. 1977) is usually two or three orders 
of magnitude higher than the observed values, while mod- 
els where thermal conduction is ignored predict a Lx that is 
typically two orders of magnitude lower than observed. Most 
of the time the descriptions are rather simple, they consider 
for instance that the cluster wind interacts with an homoge- 
neous and isotropic ISM, and neglect the effects of radiative 
cooling. For M17 this approach is not adequate because it 
is a star forming region surrounded by asymmetrically dis- 
tributed high density molecular clouds. 

In this work we present three-dimensional numerical 
simulations based on the properties of M17. The simulations 
were carried out with the adaptive grid code YGUAZU-A, 
and considered the main characteristics of the ISM that sur- 
rounds M17. Synthetic thermal X-ray emission maps were 
obtained from the numerical results in order to compare the 
morphology and the total luminosity directly with ROSAT 
X-ray observations (with a field of view of ~ 2°). 

This manuscript is organised in the following way: in 
section 2 we explain the modelling of M17, the initial condi- 
tions of numerical simulations, and describes the simulation 
of the thermal X-ray emission; an archival ROSAT image of 
this superbubble is presented in section 3; the results and 
comparison with observations are given in section 4; and 
finally in section 5 we summarise our conclusions. 



2 MODELLING THE M17 NEBULA 

As it was mentioned above, the stars that belong to the 
NGC 6618 cluster are too young to generate SN explosions. 
For this reason, the M17 nebula was modelled by the inter- 
action of stellar winds alone. We only considered the stars 
that dominate the mechanical luminosity of the cluster, they 
have high terminal velocities and mass loss ra tes, the values 
are ba sed on the spectral types reported by lHanson et al.l 
(1997), we list them in Table [T] The values for the terminal 
velocity and mass loss rate correspond to the upper limits 
reported in the literature (de Jager et al. 1988; Wilson & 
Dopita 1985; Leitherer 1988; Prinja et al. 1990; Lamers & 
Leitherer 1993; Fullerton et al. 2006). This has been done 
with the intention to increase the X-ray luminosities ob- 
tained from our simulations, and to compensate for the fact 
that all the other wind sourc es have been disregarded. Fur- 
thermore, a recent study of Hoffmcis ter et alT (|2008l ) has 
shown that the number of stars in the NGC 6618 cluster 
is even higher than previous estimates. 

For the 3D numerical simulations we consider that the 
plane of the sky corresponds to the xy— plane of our simu- 
lation. Table [T] gives the position of the stars in equatorial 
coordinates (J2000), which can be translated to parsecs con- 
sidering that the cluster is at a distance of 1.6 kpc. 

Since we do not know the individual line-of-sight dis- 
tance (z— coordinate) to the stars, we produced three dif- 
ferent realisations of randomly picked positions in z, while 
keeping the same xy configuration (see Figure The max- 
imum of the distribution from which the z positions were 
sampled was set to the maximum separation in the plane of 
the sky, and the mean was set to the distance to the clus- 
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Table 1. OB Stars in M17 and their stellar wind parameters. 



Name 


R.A.(J2000.0) 


DEC.(J2000.0) 


Optical 


K-Band 




M 






m 


C ' ") 


Spectral type 


Spectral type 


( km s" 1 ) 


CM© yr- 1 ) 


References 


B98 


18 20 35.45 


-16 10 48.9 


09 V 


k09-Bl 


2200 


2.6 x 10" 7 


1,2,3 


Bill 


18 20 34.55 


-16 10 12.1 


05 V 


k03-04 


3250 


1.86 x 10" 6 


3 


B137 


18 20 33.14 


-16 11 21.6 




k03-04 


3370 


2.4 x 10" 6 


4,5 


B164 


18 20 30.92 


-16 10 08.0 


07-08 V 


k07-08 


3015 


2 x 10~ 6 


1,4 


B174 


18 20 30.54 


-16 10 53.3 




k03-O6 


3370 


6.5 x 10~ 7 


1,4 


B181 


18 20 30.30 


-16 10 35.2 




k09-B2 


2200 


2.6 x 10~ 7 


1,2,3 


B189 


18 20 29.92 


-16 10 45.5 


05 V 


k03-04 


3250 


1.86 x 10" 7 


3 


B260 


18 20 25.94 


-16 08 32.3 


07-08 V 


k03-04 


3015 


2 x lO" 6 


1,4 


B289 


18 20 24.45 


-16 08 43.3 


09.5 V 




1500 


2.5 x 10~ 7 


1,6 


B311 


18 20 22.76 


-16 08 34.3 




k09-B2 


2200 


2.6 x 10~ 7 


1,2,3 


OI 345 


18 20 27.52 


-16 13 31.8 


06 V 


k05-06 


3065 


6.5 x 10~ 7 


1,4 



References.-(l) de Jager et al. 1988; (2) Wilson & Dopita 1985; (3) Leitherer 1988; (4) Prinja et al. 1990; 
(5) Lamers & Leitherer 1993; (6) Fullerton et al. 2006 



ter (1.6 kpc.) The resulting yz distributions are shown in 
Figured 

The NGC 6618 cluster is bounded by two molecular 
clouds at the North and South- West, observed in HI (Felli 
et al. 1984). These clouds have the appearance of a wedge 
which confine the resulting stellar cluster wind, and produce 
an elongated structure to the East. We model the clouds as 
two bars one of them is horizontal (we will reefer to it as 
the 'top bar') and the other one tilted 45°in the xy-plane 
(which we will call the 'bottom bar'). Their surfaces are flat, 
and cover the entire extent of the computational domain 
along the z axis (i.e. the geometry displayed in Figure [2] is 
the same for all values of z). The bars have a high density 
contrast with respect to the surrounding ISM (of two orders 
of magnitude). The stars (wind-sources) were placed inside 
the wedge formed by the two bars, as shown schematically 
in Figure [2] 



2.1 Initial setup 

We have carried out the 3D nume rical simulations with th e 
YGUAZU-A hydrodynamical code (Ra ga et al ] l2000l . l2002h . 
The code integrates the gasdynamic equations with a 
second-order accurate implementation of the flux vector 
splitting method of van Leer (1982) on a binary adap- 
tive grid. This code was first tested comparing simula- 
tions of shock waves and laser laboratory experiments 
(one bubble shock wave, Raga et al. 2000; two inter- 
acting bubbles, Velazquez et al. 2001). In last 8 years, 
the YGUAZU-A code has been widely used to model sev- 
eral diffe rent astrophysical scenarios, such as Her big Haro 
objects (|Raga et al.l |2002|; iMasciadri et alj |2002|). super- 
novae remnants l)Schneiter et al.l l2006: Rcvcs-Iturbidc ~et al.l 
2008), .jets in plane tary nebulae (jVelazauez et al.l 120041; 
Guerrero et al. 20081). interactin g win ds (jGonzalez et al.l 
20041; Rodrigue z-Gonzalez et all 120071 ). and exoplanets 
l|Schneiter et al.ll2007l )~ 

For this problem we use a computational domain with 
a physical size of (12 x 12 x 4) pc (as— ,y— , and z— axis, 
respectively), and allow five levels of refinement, such that 
the resolution at the finest level of is 7.27 x 10 16 cm. This 
would correspond to 512 x 512 x 128 pixels in a uniform grid. 




Figure 2. Stellar positions in the plane sky (the xy— plane). 
These positions were calculated considering the R.A and Dec. 
positions listed on Table [T] and a cluster distance of 1.6kpc. The 
shaded area represents the location of the dense molecular clouds 
that surround the M NGC 6618 cluster. 



We chose this spatial resolution (3" if a distance of .1.6 kpc 
to M17 is assumed) to match the resolution of ROSAT raw 
observations (~ 4", see §3.) 

The stellar winds are imposed in spheres (centred at 
the stellar positions) of radius R w = 4.4 x 10 17 cm, which 
corresponds to 6 pixels at the maximum resolution of the 
adaptive grid (always present at the wind sources). Within 
these spheres, we inject at every time-step material at 
T w — 1000 K, with the outward velocity given in Table [T] 
for each star. The density inside the spheres follows an r~ 2 
law (where r is the radial coordinate measured outwards 
from the stellar position), scaled to yield the mass loss rates 
in Table □ 

The stars and the bars are located in the top-right quad- 
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j v (n,T) oc n 2 . To produce X-ray maps the absorption is 
included, assuming a uniform hydrogen column density of 
N H = 3.2 x 10 21 cm -2 . 



3 ROSAT IMAGE OF M17 

M17 was observed with the ROSAT satellite, with the Posi- 
tion Sensitive Proportional Counter (PSPC) detector. This 
instrument is sensitive to X-ray photons with energies in the 
range 0.1 — 2.4 keV and has an spectral resolution energy of 
~ 40% at 1 keV, with a field of view of ~ 2°, covering the 
entire dimensions of superbubble M17 (with an angular size 
of ~ 23' x 20'). 

We have used an image (binned to 4" per pixel) 
from archival ROSAT Catalogue of PSPC WGA sources, 
which corresponds to the PSPC observations of M17 
(id:WP500311) obtained on September 12-13, 1993, with 
an exposure time of 6.7 ks. In order to study the spatial 
distribution of the X-ray emission from M17, we have in- 
creased the signal to noise ratio using the standard proce- 
dure of smoothing the image (Dunne et al. 2003, Townsley 
et al. 2003). To do this, we convolved the observations with 
a Gaussian function with an effective PSF of 5.4 pixels. In 
Figure [5] we show the original (raw data) image and the 
smoothed version. 
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Figure 3. Stellar configuration in the yz— plane for each of the 
models, as indicated in the legend. 



rant of the numerical domain, in the xy— plane. The bars 
were initialised with a uniform number density of 1000 cm -3 
and a temperature of 400 K in order to stay in pressure equi- 
librium with an ISM that has a density of 10 cm -3 , and a 
temperature of 10 4 K. 



2.2 Simulated X-ray emission 

The numerical simulation provides us with density and tem- 
perature distributions for each model, which we combine 
with synthetic X-ray spectra to simulate the X-ray emission 
of M17. 

The s ynthetic spectra were obtained wit h the chianti 
database (jDere et al.ll 19971 ; lLandi et al.ll2006h . We calculate 
the X-ray emission coefficient j v (n, T) in the energy range of 
0.1-2.4 keV to coi ncide with that of ROSAT. Based on ob- 
servational works (Dunne et al.ll2003l ; iTownslev et al.ll2003l ; 



iGarcia-Roias et~aL 2007 ) a solar abundance was assumed 
fo r MPij], togethe r with the ionisation equilibrium model 
of iMazzotta et all fl998). The j u (n,T) coefficients are cal- 
culated in the limit of low density regime, which results in 



1 In §4 section we discuss the effect of the metalicity on the total 
X-ray luminosity. 



4 RESULTS AND DISCUSSION 

We ran the models up to an integration time of 1.3 x 10° yr, 
which is the estimated superbubble age (Weaver et al. 1977). 
In Figure [4] we present X-ray emission maps for the three 
models at their latest integration time. These maps show 
the absorbed X-ray emission, assuming a column density of 
3.2 x 10 21 cm -2 as found from ROSAT observations (Dunne 
et al. 2003). 

The interaction of many individual stellar winds 
will coalesce into a common cluster wind , as st udied 
for massive stellar clust er s by ICanto et al.l (|2000t ) and 
iRodrfguez-Gonzalez et al] l|2007h . However, unlike those 
models that consider a uniform ISM, we can see from Fig- 
ure |3] that the resulting X-ray emission is not homogeneous. 
All models show a strong concentration of X-ray emission in 
regions to the top right of each map, especially for model 2. 
This enhancement corresponds to strong shock waves pro- 
duced by the collision of the cluster wind with the dense 
molecular cloud or bars. These bars are flat and focus the 
cluster wind flow and produce the elongated morphology ob- 
served in X-rays. To the opposite side the cluster wind can 
expand more freely producing through a 'champagne'-type 
flow a superbubble of ~ 12 pc in size. From Figure 3] the 
structure seen in our simulations at the position x ~ lOpc, 
y ~ 7pc corresponds to the wind of an isolated star embed- 
ded in the the southern bar, which gets confined due to the 
high density of its surrounding medium (see Fig. [2}. 

The surrounding dense bars have another important ef- 
fect, namely to increase the total X-ray luminosity Lx- This 
effect depend on the geometry of the bars, which in our 
models we have considered to be flat (independent of z). If 
they would have a concave section, the cluster wind would 
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Figure 4. Simulated X-ray emission maps obtained from models 
1, 2, and 3 (from the top to the bottom) in the energy range 
0.1 - 2.4 keV, an absorption due to N H = 3.2 X 10 21 cm~ 2 has 
been included. Maps corresponding to models 1 and 2 show more 
evolved super bubbles (produced by the cluster wind). The bar on 
the top gives the logarithmic grey-scale in units of erg cm -2 s — 1 
sr _1 . The vertical and horizontal axis are given in units of pc. 

be more strongly focused, increasing the Lx, and the oppo- 
site would be true if they were convex. To analyse influence 
of the bars, we have run an additional model (not presented 
here), removing them from the numerical setup. In that case 
the cluster wind is free to sweep the less dense ISM, and re- 
sults in a Lx 2 orders of magnitude lower than the obtained 
one from Models 1, 2, and 3. This shows the importance of 
considering the characteristics of the environment of M17, 



in particular the dense molecular cloud that surrounds the 
open cluster NGC 6618. 

In Table [2] we list the Lx obtained for all the mod- 
els. We must notice that the X-ray luminosity listed there 
has been calculated directly from the X-ray emission coef- 
ficients derived from the simulations, thus no correction for 
absorption is needed. In other words, it should be compared 
with the X-ray luminosity of observations after these obser- 
vations have been corrected for absorption. Models 1, 2, and 
3 produce similar luminosities, all of them on the order of 
60 per cent of the observed value (although Model 1 gives 
a slightly lower value). One has to bear in mind that the 
gas in the bar is not uniform in reality, but it has high den- 
sity clumps (Townsley et al. 2003). The inclusion of such 
clumps could somewhat increase the total X-ray luminosity 
compared with that obtained in our models. 

Other factor that can change the obtained Lx is the 
metal abundance. To estimate importance of this factor, we 
generated new X-ray coefficients (with the CHIANTI atomic 
database) considering different values of the metalicity Z 
around of the solar value. The resulting luminosities are 
listed in table For Z = 0.3, the X-ray luminosity is 2.7 
times lower than the obtained for Z — 1 (solar abundance). 
We must note, however, that this is only a rough estimate. 
A detailed study of the effect of metalicity on the X-ray 
luminosity would have to include the dynamical effect that 
result s from a modified radiativ e cooling. In any case, recent 
work ( Garci'a- Roias et al1l2007l ) show that abundances close 
to the solar ones are adequate for M17. 

Dunne et al. (2003) obtained a total X-ray luminosity 
of 5 x 10 35 erg s _1 for M17, employing both the Weaver et 
al. (1977) bubble model, and the Chu et al. (1995) super- 
bubble model. Such value is two orders of magnitude higher 
than the observed one. These models are based on the evo- 
lution of bubbles (and superbubbles) into an homogeneous 
medium, and they are adiabatic, although they include ther- 
mal conduction effects. In our case, it was not necessary to 
include the last physical process because we have obtained 
a Lx already compatible with observations, taking only into 
account the inhomogeneity of surrounding ISM, and includ- 
ing the radiative cooling. 

Although thermal conduction effects should be present, 
as firstly introduced in the wind-blown bubble model by 
Weaver et al. (1977), there are important discrepancies be- 
tween the X-ray luminosities predicted by that model and 
what is observed. For instance, X-ray emission should be 
detected in almost all the known wind-blown bubbles and 
superbubbles, but this is not the case (see Mac Low 2000, for 
a discussion). In particular, for M17, a direct application of 
Weaver's et al. model (which considers thermal conduction 
as a crucial ingredient for the X-ray emission of the hot in- 
terior) predicts an X-ray luminosity about 100 times higher 
than what is measured (Dunne et al. 2003) . This leads us to 
think that perhaps some mechanisms are present that inhibit 
thermal conduction, at least in the form that it is treated 
in Weaver's et al. model. The presence of magnetic fields is 
often invoked as one of the possible inhibiting factors. In our 
particular case, we show that whatever mechanism might be 
at work, thermal conduction do not seem to play an impor- 
tant role in the prediction of X-ray emission from the M17 
superbubble. 

The X-ray emission distribution for all models (see Fig- 
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Table 2. Total X-ray luminosities. 



L x (10 33 erg s" 1 ) 
ROSAT observation^ 2.5 
Model 1 1.4 
Model 2 1.6 
Model 3 1.6 



a Dunne et al. 2003 

Table 3. Mctalicity vs Total X-ray luminosities. 



Z L x (10 33 erg s' 1 ) 

0.3 0.6 

0.5 0.9 

0.7 1.2 

0.9 1.4 

1.0 1.6 

1.1 1.7 
1.3 2.0 



ure[4]) show a general gradient, with lower emission to the 
left, in coincidence with the observations (see Figure [5}. Also 
in our models, on the left of the X-ray emitting material, 
strong emission is observed in form of filaments (see Fig- 
ure These filaments trace the shock front of champagne 
flow described above, propagating into the surrounding ISM. 
In ROSAT observations, bright regions or 'clumps' can be 
seen to the East (corresponding to the left, in our simula- 
tions) . These clumps can be produced by the propagation of 
strong shock waves into an inhomogeneous medium. How- 
ever, these features are smoother than those obtained from 
the simulations probably due to the difference between ob- 
served and simulated resolutions. To make a fair comparison 
we show in Figure [6] a smoothed synthetic map of model 2, 
produced by convolving the numerical results with a Gaus- 
sian beam in order to reproduce the same resolution of the 
ROSAT (smoothed) image in Fig. [5] 



5 CONCLUSIONS 

We carried out 3D numerical simulations using the adaptive 
grid code YGUAZU-A to model the M17 superbubble. Ra- 
diative losses have been included in the simulations. Three 
different runs of the same model were made using the known 
positions of the dominant wind sources on the plane of the 
sky (corresponding to the xy positions in our Cartesian 
grid), and considering three different distributions along the 
line of sight (aligned with our 2— axis). 

Our results show that the inclusion of the main fea- 
tures of the surrounding ISM (i.e. the presence of two dense 
bars or clouds around the M17 stars) plays a crucial role to 
explain both the observed morphology and the total X-ray 
luminosity of this object. 

On the one hand, the bars produces in all models a 
'champagne' flow effect, in which the resulting stellar cluster 
wind is focused by such bars, producing an elongated shape 
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Figure 5. ROSAT image of superbubble M17, obtained from 
archival data (ROSAT catalogue of PSPC WGA sources) in the 
energy range 0.1 — 2.4 keV. Top panel: raw data binned to 4" per 
pixel. Bottom panel: Smoothed image obtained after convolving 
with a Gaussian function with a PSF of 5.4 pixels. 

in the x— direction, thus reproducing the elongated shape 
in the East- West direction of the observations. 

On the other hand, the collision of the cluster wind 
with the bars maintains a high-temperature, high-density 
distribution in the region between the bars. This is a product 
of the multiple reflected shocks, and enhances the total X- 
ray luminosity (Lx). We obtain values of the total X-ray 
luminosity of ~ 1.6 x 10 33 erg s _1 , which are in agreement 
with observations. 

The different z— stellar distributions employed in our 
runs had a rather small effect on both the morphology and 
Lx- This could be somewhat expected, because the same 
mechanical luminosity input is used for all the runs, and 
because away of the cluster centre the common wind that 
forms should be similar. Finally, it is noticeable that we 
obtain a reasonable estimate of the total X-ray luminosity 
without the need of including thermal conduction effects. 
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